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Abstract 

We study a problem of failure of two interdependent networks in the case of correlated degrees 
of mutually dependent nodes. We assume that both networks (A and B) have the same number of 
nodes N connected by the bidirectional dependency links establishing a one-to-one correspondence 
between the nodes of the two networks in a such a way that the mutually dependent nodes have 
the same number of connectivity links, i.e. their degrees coincide. This implies that both networks 
have the same degree distribution P(k). We call such networks correspondent^ coupled networks 
(CCN). We assume that the nodes in each network are randomly connected. We define the mutu- 
ally connected clusters and the mutual giant component as in earlier works on randomly coupled 
interdependent networks and assume that only the nodes which belong to the mutual giant com- 
ponent remain functional. We assume that initially al—p fraction of nodes are randomly removed 
due to an attack or failure and find analytically, for an arbitrary P(k), the fraction of nodes fi(p) 
which belong to the mutual giant component. We find that the system undergoes a percolation 
transition at certain fraction p = p c which is always smaller than the p c for randomly coupled 
networks with the same P(k). We also find that the system undergoes a first order transition at 
p c > if P(k) has a finite second moment. For the case of scale free networks with 2 < A < 3, the 
transition becomes a second order transition. Moreover, if A < 3 we find p c = as in percolation 
of a single network. For A = 3 we find an exact analytical expression for p c > 0. Finally, we find 
that the robustness of CCN increases with the broadness of their degree distribution. 
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I. INTRODUCTION 



The robustness of interdependent networks has been recently studied by Buldyrev et al. 
[H within the framework of the mutual percolation model. They found that two randomly 
connected networks with arbitrary degree distributions randomly coupled by bidirectional 
dependency links completely disintegrate via a cascade of failures if the fraction p of the 
nodes which survive the initial attack is less than some critical value p = p c > 0. Moreover, 
the transition at p c is of the first order type, i.e. the fraction of the functional nodes p,(p) 
which survive after the cascade of failures has a step discontinuity at p = p c changing from 
A*c — MPc) > f° r P = Pc to zero for p < p c . This behavior was observed even for scale 
free (SF) networks with a power law degree distribution P(k) ~ k~ x with 2 < A < 3. 
The explanation of this behavior is based on the fact that in that model the nodes with 
large degree (hubs) in one network may depend on the nodes with small degree in another 
network. The nodes with small degree can be isolated from a giant component in one 
network by removal of a small fraction of nodes and thus cause the malfunction of the 
hubs in the other network. In real world interacting networks, the hubs in one network are 
more likely to depend on the hubs of another networks. This can significantly enhance the 
robustness of the interdependent networks. In general, the correlations among the degrees 
of the mutually dependent nodes can be described by a matrix P(ki\k2) which specifies the 
conditional probabilities to find a node with degree k\ in one network, provided it depends 
on a node with degree k^ in another network. This matrix can be quite complex and may 
depend on many parameters. For each parameter set the model can be readily studied by 
computer simulations, but in order to get general understanding of the correlation effects it 
is desirable to solve the problem analytically at least in some limiting cases. 

In this paper we study the mutual percolation problem in the case of the strongest pos- 
sible correlations, namely we studied the case in which both networks (A and B) have the 
same number of nodes N connected by bidirectional dependency links establishing a one- 
to-one correspondence between the nodes of the two networks in such a way that mutually 
dependent nodes have the same number of connectivity links, i.e. their degrees coincide, i.e. 
P{k\\k2) = 1 for k\ = &2 and P{k\\k2) = otherwise. This implies that both networks have 
the same degree distribution P(k). For brevity we will call such networks correspondent ly 
coupled networks (CCN), while we will refer to the model studied in Ref. [ij as randomly 
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coupled networks (RCN). As in Ref. we assume that the nodes in each network are 
randomly connected. We define the mutually connected clusters and the mutual giant com- 
ponent as in Ref. [l| and assume that only the nodes which belong to the mutual giant 
component remain functional. 

We assume that initially a 1 — p fraction of nodes are randomly removed due to an attack 
or failure and find analytically the fraction of nodes p,(p) which belong to the mutual giant 
component. We find that as in Ref. [l], the system undergoes a percolation transition at 
certain fraction p = p c which, however, is always smaller than the pa for RCN with the same 
degree distribution with the exception of random regular graphs 2] for which both values 
coincide. Moreover, we find that the system undergoes a first order transition at p c > 
if the degree distribution has a finite second moment. For the practically imqportant case 
of SF networks {3]- 7] with 2 < A < 3, for which the second moment diverges, the transition 
becomes a second-order transition. If A < 3 we find that p c = as in percolation of a single 
network [8], while for A = 3 we find an exact analytical expression for p c > 0. The change 
in transition order has been also observed in interdependent networks with partial coupling 
9J. We also investigate how the broadness of the degree distribution affects the robustness 
of CCN. 

II. GENERATING FUNCTIONS AND THE CASCADE PROCESS 
A. First stage 

We will describe the stages of the cascade of failures in CCN in terms of the generating 



function of their degree distribution [10|, 111 ] 



G(x) = J2P(k)x k (1) 

k=0 



and the generating function of the associated branching process 12] 



{X) ~ G'(l) ~ (k) dx ' { ) 

where (k) = G'(l) is the average degree. It is known that the degree distribution P(k,p) 
of a network, from which a fraction of nodes 1 — p is randomly removed is related to the 
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original distribution P(k) through a binomial expansion 10] : 



p(k',p) = J2P(k)p k '(i-p) k - k 'c l 



k > 



(3) 



k>k' 



where C| = k\/\k'\(k — A;')!] are binomial coefficients. Accordingly 10], the generating 
function of this distribution is 



G(x,p) = G(xp + 1 — p). 



(4) 



The fraction of nodes which do not belong to the giant component of a network is given 



by 
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r = G(f), 

where / is a smallest nonnegative root a transcendental equation 

f = H(f). 



(5) 



(6) 



The degree distribution of nodes which do not belong to the giant component is given by 

Q 

P (kJ) = P(k)f k /r. (7) 
Accordingly the degree distribution of nodes in the giant component is given by 



P i (k,f)=P(k)(l-f k )/(l-r). 



(8) 



Thus the degree distribution in the giant component of a decimated network after random 
removal of a 1 — p fraction of nodes is 



P(k', f,p) = P(k',p)[l - f{pf']/[l - r(p)}. 



where 



r (p) = G(f(p),p) 
and f(p) satisfies the transcendental equation 

f(p) = H(f(p),p). 



(9) 



(10) 



(11) 



In order to find the original degree distribution in the giant component of network A we 
must restore the links which lead to the randomly removed nodes. If a node in the decimated 



network A has a degree kl it might have any degree k > k! in the original network A with 
probability P(k\k') given by the Bayes' formula: 

P{k\k>) = P(k)C k k 'p k '(l-p) k - k '/P(k,p). (12) 

Thus the total probability that a node in the giant component has a degree k is 

P 1 (k) = £ P{k)CW{l-p)*-« P fl f 'd (13) 
k><k P l fc >£V 

or using Eq. ([9]) 

The generating function of this degree distribution is 

G(x) - g(gfi) 
= l-G( tl ) ' (15) 

where £i = /ip + 1 — p and /i = /(p). The fraction of nodes in the giant component of the 
decimated network A is 1 — rj., where r x = G(ti), Because the decimated network has Np 
nodes, the size of the giant component A x of network A after random removal of (1 — p) 
nodes is Ni = Np(l — ri). 



B. Second stage 

We assume that only nodes which belong to Ai are functional; thus after the first stage of 
the cascades of failures, only p(l — ri) < p fraction of nodes in network B remain functional. 
Thus we expect further disintegration of network B at the second stage of the cascade and its 
giant component B2 will be even smaller than Ai. We define a set of nodes Bi = D(Ai) by 
projecting Ai onto network B using the one-to-one correspondence D between the nodes of 
networks A and B established by dependency links. Since the degree of each node in network 
B is the same as the degree of its dependent node in network A, the giant component A x 
of network A obtained at the first stage of the cascade has the same degree distribution 
as the set Bi in network B. Moreover, from the point of view of network B the nodes in 
Bi are randomly selected and randomly connected. Thus, to compute B2 we can use the 
same approach used at the first stage, but applied to the new network Bi with the new 
degree distribution given by Eq. ()14p . The only problem is that many of the links outgoing 

5 



from network Bi are ending at the nodes which do not belong to network Bi and thus 
for computation of B2 these links must be removed. The probability p\ of a random link 
originating in network B x to end up in Bi is equal to the ratio of the number of links 
originating in network B]_: 

L x = N x £ kPi(k) = pN(k)(l - G'it^h/ik)) (16) 

to the total number of links N(k). Therefore, 

Pi = ]^)=P(l-«i), (17) 

where 

s 1 = t 1 G'(t 1 )/(k). (18) 
Accordingly, the degree distribution of links connecting the nodes of network B x is 

PiW,p) = E fi(*)pf (i -Pi) k ~ k 'cl' (19) 

k>k' 

and the generating function of this distribution is 

G 1 (x, Pl ) = [G(x Pl + I - Pl ) - G^xp, + 1 - Pl ))]/{1 - n). (20) 

Thus the size N2 of the giant component B2 is N2 = p(l — ri)(l—r2)N, where r2 = Gi(f2,pi) 
and f2 = Hi(f2,Pi)- Introducing a new notation 

£2 = f?Pi + I-P1 (21) 
and taking into account Eq. fl20l we see that 

h (*)(!-„) (22) 

and N 2 = p(l - r x )[l - (G(t 2 ) - G(t 2 tx))/(1 - n)]. Using that n = G(ti), we get 

iV 2 = p[l - Gfa) - G(t 2 ) + Gihh^N (23) 

We can compute the original degree distribution P 2 (k) in B 2 using the Bayes' formula the 
same way as we obtained the distribution Pi (A;) 



C. Third stage 



On the third stage of the cascade we will compute the giant component A 3 of network 
A which is the result of further disintegration of A 1 because the nodes in A x which do not 
belong to B 2 failed at the second stage. We can again apply the same technique, with 
the only difference that now A 2 = -D(B 2 ) is not a random subset of nodes of A but they 
are taken out of its connected giant component A x . Accordingly we must find an effective 
degree distribution and the effective size of a network which would reproduce A 2 by random 
selection of nodes out of this original network. Since the degree distribution of nodes in A 1 
is 

^-^r^h < 25 > 

and the degree distribution in A 2 is 

^ H'-S-ff - (26) 

the selection of A 2 out of Ai has the same effect as random selection of a fraction p(l — r\) 
out of the entire network with degree distribution P^ik) = P{k){l — t 2 )/(l — r 2 ). Now we 
can see that the problem of the third stage is completely equivalent to the second stage with 
t\ replaced by t 2 , fi replaced by / 2 , T\ replaced by r 2 , t 2 replaced by £3, / 2 replaced by f^, 
and r 2 replaced by r 3 . 

D. Recursive relations 

Generalizing, for stage i we arrive to a recursive relation between ti and U+i. Namely, 
once we know ti we can find U+i, as well as the size of the giant component at the stage 
i + 1 

N i+1 = p[l - G(U) - G(t i+1 ) + G(tit i+1 )]N (27) 
and the degree distribution of the nodes inside this giant component 

^< fc > - p W i- G (g:gr + ^,) (28) 

In order to find t i+1 from ti we repeat the steps used deriving t 2 from t 1 by first introducing 

s i = t i G'(t i )/(k), (29) 
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and 

Pi = p(l - Si ) (30) 
in analogy to Eqs. ( 129]) and f fT7|) . Then 

= /i+lPi + 1 ~Pu ( 31 ) 
where satisfies a transcendental equation analogous to Eq. (l22l) 

fi+1 " <*>(i-*) • (32) 

Excluding and Sj from Eq. fl32|) we find that t i+1 is given by the smallest non-negative 
root of the equation 

t l+l = (l-p) + jj-r[G' (ti)ti + G'{t i+1 ) - UG'(t i t i+1 )]. (33) 

To start the iterative process we must take into account the definition of t\ = pfi + 1 — p 
and Eq. (ITT]) which is equivalent to a transcendental equation 

tl = (l-p) + JLG\t 1 ), (34) 
which is the same as Eq. ( 133]) if we introduce to = 0. 

III. THE MUTUAL GIANT COMPONENT AND THE PHASE TRANSITION 

The cascade of failures will stop when = t{ — t and hence the fraction of nodes in 
the mutual giant component p, = lim^oo Ni/N, is given by the simplified equation (127]) : 

p = p[l - 2G(t) + G(t% (35) 

where t is the smallest non-negative root of the equation 

t = (1 -p) + -^[(1 + t)G\t) - tG\t 2 )} = 1 -p[l - (1 + t)H(t) + tH(t 2 )]. (36) 
{k} 

The right hand side of Eq. ( 136]) has zero derivative at t = 1, if G"(l) is finite. This 
condition is equivalent to the existence of the second moment of the degree distribution. 
Thus one can see [Fig. [JJ[a,b)] that for finite second moment and small enough p, Eq. ( 136]) 
has only the trivial solution t — 1 corresponding to p = and, therefore, to the complete 
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disintegration of the networks. As p increases, a nontrivial solution p > will emerge at 
p = p c , at which point the right hand side of Eq. fl36|) will touch the straight line representing 
the left hand side at t — t c ; at that point the slope of both lines is equal to 1. Since at 
t = 1 the slope of the right hand side is zero, t c must be smaller than 1 and thus the 
mutual percolation transition is of the first order, where p changes form zero (for p < p c ) to 
P > Pc > (for p > p c ). The value of p c is given by Eq. ( |35|) computed at t = t c . 
An efficient way of finding p c is to solve Eq. ( 136]) with respect to 1/p: 

l-(l + t)H(t)+tH(t 2 ) = 1 

1 -t p 1 ' 

and find the maximum of the left hand side with respect to t [Fig. [2]. The left hand side of 
Eq. ( l3Tj) is a curve which changes from 1 — H (0) = 1 — -P(l)/ (k) at t = to zero at t = 1. At 
t = it has a positive slope 1 — (-P(l) + 2P(2))/ (A;)), so it must have an absolute maximum 
at t c G (0, 1). The equation for t c can be readily obtained by differentiation of Eq. (1371) 

1 - 2H(t c ) + H(t 2 c ) - (1 - t 2 c )H'(t c ) + 2t 2 c (l - t + c)H\t 2 c ) = 0. (38) 

The value of the left hand side of Eq. (1371) at t — t c gives l/p c - If the value of this 
maximum is less than 1, than the networks do not have a mutual giant component at any p. 



IV. SPECIAL CASES 

Figure [H shows the graphical solutions of Eq. (1361) for several special cases of degree 
distributions of CCN. 



A. Erdos-Renyi networks 

For Erdos-Renyi (ER) networks [2I. [15I H(t) = exp[(£;)(i — 1)] and the maximal value of 
the left hand side of Eq. (1371) monotonously increases with (k) . This can be readily seen by 
differentiating Eq. fl37j) with respect to (k). The maximal value reaches 1 at (k) = 1.706526, 
below which correspondently coupled ER networks disintegrate even without any initial 



attack or failure (" 
networks is 2.4554 



? ig. [2]). Note that the equivalent value of (k) for randomly coupled ER 
II- 
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B. Random regular graphs 



For RR graph [Fig. Hfb)] in which all the nodes have the same degree k = (k), G(t) = 
and H(t) = t<*> -1 . Then t satisfies 

t = (l-p)+ p{t^- 1 + t {k) - t 2 ^" 1 ), (39) 

and 

fi = p{\ - t<*>) 2 . (40) 

Equations (1391) and (140]) can be obtained by simpler methods presented in Ref. jlj for RCN, 
since for the case of random regular graphs, the degrees of all the nodes in both networks 
coincide and therefore the CCN and RCN models are equivalent. Indeed from Eq. (1) of 
Ref. [l| it follows in a special case of coinciding degree distributions of the coupled networks 
that 

V = p(l-G(t)) 2 , (41) 

where 

t = l-p(l-G(t))(l-H{t)). (42) 

If G(t) = t (fc) and H{t) = t^' 1 , Eqs. gTJ and O are equivalent to Eqs. flUD and 
respectively. 



C. Scale free networks 

For scale free (SF) networks with A < 3 [Fig. 00(c)], the second derivative of the right 
hand side of Eq. (1361) is infinite at t = 1, which means that a nontrivial solution exists at 
any p > since in the vicinity of t = 1 the straight line representing the left hand side of Eq. 
fl36|) is always above the curve representing the right hand side, while for t — 0, the curve 
is always above the line. This means that SF CCN are as robust as a single SF network for 
which p c is always zero. 

For the marginal case of A = 3 [Fig. [Tj(d)] G"(t) diverges as ln(l — t) when t — > 1 and thus 
the left hand side of Eq. (I3"B]) has a finite slope at t — 1. Accordingly p = p c > but the 
nontrivial solution emerges at t c = 1, so the transition becomes of the second order. For the 
case of P(k) = [k min /k) 2 - {k min /(k + l)) 2 for k > k min = 1,2, .. and P(k) = for k < k min 
we can find p c analytically. Indeed, in this case P(k) behaves asymptotically as 2fc 2 rain /fc 3 . 
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For k — y oo the leading term in G"(t) becomes 1k 2 min t k jk\ so, G"(t) = —2k 2 min ln(l — t) + c(t), 
where c(t) is a continuous function. Accordingly, the slope of the right hand side of Eq. fl36|) 
at t — 1 becomes p4/c^ n ln(2)/(/c), where 



The critical threshold is thus 



1 I 7T 2 V^fcmm _1_ 

Pc " 41n(2) 1 ' 



For A; min = 1 p c = 0.59328456 and for k min = 2 p c = 0.32277924. 



D. Effect of the broadness of the degree distribution 

It follows from Fig. [1] that for the same (k) = 3, p c of the RR networks (0.758751) is 
greater than the p c of the ER networks (0.6499451). Moreover, for SF networks with A = 3 
and k m i n = 1, for which the average degree is ir 2 /6 < 3 we have even smaller p c = 0.59328456. 
For SF networks with A = 3, (k) — 3 we can estimate p c = 0.35 which is much smaller than 
the p c for the narrower ER and RR degree distribution. For SF networks with A < 3, which 
are even broader, p c = for any (A;). This is in a complete agreement with the trend observed 
in percolation of single networks, for which the robustness increases with the broadness of 
the degree distribution if one keeps (k) constant, but is opposite to the trend observed in 
Ref. |l| for RCN. 

In order to investigate the effect thoroughly, we study several classes of degree distri- 
butions for a number of values of (k). Figure [3] shows p c as function of (k) for RR, ER, 
uniform, and SF with A = 3 degree distributions. For each value of (k) the variance of 
SF degree distribution (oo) is greater than the variance of the uniform degree distribu- 
tion ((k)({k) + l)/3), which is greater than the variance of ER degree distribution (k) 
which is greater than the variance of RR degree distribution (0). Indeed, Fig. [3] shows 
that Pc(SF) < p c (uniform) < p c (ER) < p c (RR). Thus our numerical results suggest that 
CCN become more robust if their degree distribution becomes broader (provided the average 
degree is constant). This behavior is the opposite of the behavior of RCN. 

However, in general, if the measure of broadness is simply the variance of the degree 
distribution, our statement is incorrect. It is possible to find two distributions with the same 
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variances and average degrees, which have different values of p c . One particular example is 
the following two distributions Pi(0) = 0, P X {1) = Pi (2) = Pi (3) = Pi (4) = Pi(5) = 1/5 and 
P 2 (0) = P 2 (3) = 0,P 2 (1) = P 2 (5) = 1/6,P 2 (2) = P 2 (4) = 1/3 which have p c respectively 
0.683099 and 0.683657. 

V. GENERAL IMPLICATIONS ON THE NETWORK ROBUSTNESS 

Finally, we will compare the robustness of CCN and RCN with the same degree distri- 
butions. We will show that (i) the value of p c for CCN is always (except for RR networks) 
smaller than the p c for RCN and (ii) for the same p, the value of the mutual giant component 
for CCN is always (except for RR networks) larger than for RCN. 

Eq. fT4"2"l) for the randomly coupled networks can be rewritten as 



The critical value of p c for randomly coupled networks can be determined as the inverse 
maximal value of the left hand side of Eq. ( 145]) . Our proposition (i) is an obvious corollary 
of the following proposition (Hi): for any t £ [0, 1] the left hand side of Eq. fl37|) is greater or 
equal than the left hand side of Eq. ()45l) (Figj2j) . Subtracting Eq. ( )45l) from Eq. (J37l) and 
applying relation ([2]) between G(t) and H(t) we see that the inequality stated in proposition 
(Hi) is equivalent to 



We will prove Eq. (|46j) using mathematical induction. We see that for RR graphs for 
which the degree of every node is equal to m, i.e. P(m) = 1, Eq. PS"]) is satisfied as an 
equality. Assuming that it is satisfied for any degree distribution such that P(k) = for 
k < m and k > n > m, we will show that it is also satisfied for the degree distribution 
P(k) = (1 — b)P(k) for any k except for k — n + 1, for which P(n + 1) = b > 0. The 
generating function for this new distribution is obviously G = (1 — b)G + bt n+1 . After 
elementary algebra we can see that 



[1 - H(t)][l - G(t)} _ 1 
1-t p 



(45) 



tG'(t 2 ) - tG'(t) + G(t)G'(l) - G(t)G'(t) > 



(46) 



tG'(t 2 ) - tG'(t) + G(t)G'(l) - G(t)G'(t) = 
(1 - b)[tG'(t 2 ) - tG'(t) + G(t)G'(l) - G(t)G'(t)} + 
b(l - b)(G(t) - t n+1 )[(n + 1 - G'(l))(l - t n ) + G'(t) - G'(l)t% 



(47) 
(48) 
(49) 
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which proves Eq. (146]) for G provided it is true for G, if we take into account the obvious 
inequalities n + 1 > G'(l), 1 > t n , G(t) > t n+1 and G'(t) > G'{l)t n for any t e [0, 1]. This 
concludes the proof of propositions (Hi) and (i). Note that the equality sign in the above 
inequalities and hence in inequality ( 146]) is realized only for t = 1 and t = (if -P(O) = 0). 
Hence proposition (i) always implies rigorous inequality except for the case of RR graphs. 

To prove the proposition (ii) we first notice that the smallest positive root of Eq. (l3Tj) . t 1; 
is always smaller than the smallest positive root t 2 of Eq. (145]) . This a is direct consequence 
of proposition (Hi). Also we notice that the right hand side of Eq. (135 p is a monotonously 
decreasing function of t. This can be shown by differentiation and comparing the terms of 
G'(t) and tG'(t 2 ) corresponding to the same P(k), namely kP(k)t k ~ 1 > kP(k)t 2k ~ 1 . Thus 
fi(ti) > fi(t 2 )- Finally we state proposition (iv): for the same value of t, the right hand 
side of Eq. (135]) . fi(t), is greater or equal than the right hand side of Eq. (jUJ), /i r (t). One 
can prove this proposition using the same induction method we used to prove proposition 
(Hi). Combining these two results fi(ti) > n(t 2 ) > ^(h), which concludes the proof of 
proposition (ii). 

Thus CCN are more robust than RCN with the same P(k), but are still prone to cascade 
failures and, then, to first order disintegration (only if G"(l) < oo) as in case of randomly 
coupled networks. 

VI. SUMMARY 

In this work we have studied the problem of failure of CCN, i.e. coupled networks 
with coinciding degrees of mutually dependent nodes. We derive new recursive equations 
[Eqs. ( 133]) and ( 1271) ] describing the cascade of failures that are different from the analogous 
equations for RCN studied in Ref. We also find equations for the size of the mutual giant 
component [Eqs. (135]) and ( 136]) ]. as well as the efficient way of finding the critical fraction of 
nodes p = p c which must survive the initial random failure for the mutual giant component 
not to vanish, by finding the maximum of Eq. ( I37j) . 

We show that if the second moment of the degree distribution is finite, CCN disintegrate 
in a cascade of failures via a first order transition at which the mutual giant component 
suddenly drops from a positive fraction above p c > to zero below p c . This behavior is 
analogous to the behavior of RCN, with the only difference that RCN disintegrate via a first 
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order transition even when the second moment of their degree distribution diverges. 

Moreover, we show that CCN are always more robust than RCN with the same degree 
distribution. In particular, we show that scale free CCN with A < 3 disintegrate via a second 
order phase transition in the same way as non-interacting networks and thus are very resilient 
against random failure. Namely, the mutual giant component for these networks exists at 
any p > 0, but becomes infinitely small as p — > 0. Finally CCN become more robust if 
their degree distribution becomes broader (provided the average degree is constant). This 
behavior is the opposite of the behavior of RCN. 

All our analytical predictions are confirmed by simulations of coupled networks with large 
number of nodes (N > 10 6 ). 

Based on our findings we conjecture that coupled networks with any positively correlated 
degrees of mutually dependent nodes (and not just the present case of fully coincidental 
degrees) are more robust that their randomly coupled counterparts studied in Ref. JJ This 
can be attributed to the fact that the correlation between the degrees of nodes suppresses 
(or attenuates) the phenomenon of hubs becoming more vulnerable by being dependent on 
low degree nodes in a coupled network. 
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FIG. 1: Graphical solution of Eq. (j36H for various special cases of CCN. (a) ER networks with 
average degree (k) = 3. One can see that the black curves representing the right hand side of 
Eq. (|36p for different p have zero slopes at t = 1. The relevant solutions for t are given by the 
lower intersection points of these curves and a straight line y = t representing the left hand side 
of Eq. ()36p . For p = 1, this solution t = 0.0602 is indicated by a vertical straight line. The 
intersection of this vertical line with the plot of Eq. (|35p (red curve) gives the mutual giant 
component fi = 0.931. The critical p = p c = 0.6499451 corresponds to sudden disappearance of 
the nontrivial solution, (b) RR networks with k = 3. Note that for p = 1 the nontrivial solution is 
t = which means that [i = 1. The value of p c = 0.758751 is grater than the p c for ER networks 
with the same average degree shown in panel (a), (c) Analogous plot for scale free networks with 
A = 2.5. It shows that the slope of the curves is infinite for t —■ 1. One can see that in this case 
the nontrivial solution exists for any p > 0. However as p — > 0, the nontrivial solution t — >• 1, and, 
accordingly, fi — > indicating the second order transition at p = p c = 0. (d) The marginal case of 
A = 3. The slopes of the curves for t — > 1 are finite. This means that there is a critical p = p c > 
at which the slope of the curve becomes equal to 1 at t — > 1. For the displayed case of k m in = 1, 
Eq. (|44p yields p c = 0.59328456. The nontrivial solution smoothly approaches 1 as p — > p c . This 
again implies ji — > (second order transition). 
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FIG. 2: Graphical solution of Eq. (J37J for ER networks with different degree < k > illustrating 
the method of finding p c . The red curve corresponds to (k) = 3 studied in FigQJa). As < k > 
decreases below 1.706, the nontrivial solution corresponding to p < 1 disappears. We also show the 
behavior of the analogous equation (|45|) for (k) = 1.706 for RCN. In agreement with proposition 
(Hi) this curve is always below the curve with the same average degree for networks with coinciding 
degrees studied here. 
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<k> 

FIG. 3: The values of p c versus (k) for several degree distributions of increasing broadness, 
namely RR, ER, uniform and SF with A = 3. We define the uniform distribution as follows: 
P(k) = l/(2(k) + 1) for k = 0, 1, 2(k) and P(k) = for k > 2{k). For SF distribution we use 
Eqs. (03|) and while for other distributions we numerically solve Eq. (f38|) and use Eq. ([37]) 
to find p c . One can see that p c decreases ( and hence the robustness increases) with the broadness. 
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